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Abstract. The measurement of the pion form factor and, more generally, of the cross section for electron- 
positron annihilation into hadrons through the radiative return has become an important task for high 
luminosity colliders such as the <?>- or B-meson factories. For detailed understanding and analysis of this 
reaction, the construction of a Monte Carlo program, PHOKHARA, has been undertaken. Version 2.0 was 
based on a next-to- leading order (NLO) treatment of the corrections from initial-state radiation (ISR). In 
the present paper a further extension of PHOKHARA (version 3.0) is described, which incorporates NLO 
corrections to final-state radiation (FSR). The impact of combined ISR and FSR on various distributions 
is investigated and methods are presented, which will allow the extraction of the form factor, and even give 
access to inclusive photon emission due to FSR. The dependence of the results on the model for FSR is 
discussed and the impact of this contribution on the anomalous magnetic moment of the muon is evaluated. 



1 Introduction 

Electron-positron annihilation into hadrons is one of the 
basic reactions of particle physics, crucial for the under- 
standing of hadronic interactions. At high energies, around 
the Z resonance, the measurement of the inclusive cross 
section and its interpretation within perturbative QCD ^ 
E] give rise to one of the most precise and theoretically 
founded determinations of the strong coupling constant 
a s [3]. Also, measurements in the intermediate energy re- 
gion, between 3 GeV and 11 GeV can be used to determine 
a s and at the same time give rise to precise measurements 
of charm and bottom quark masses |U . The low energy re- 
gion is crucial for predictions of the hadronic contributions 
to a^i the anomalous magnetic moment of the muon, and 
to the running of the electromagnetic coupling from its 
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value at low energy up to Mz (for reviews see e.g. 
|7||5]|§]: the most recent experimental result for a M is pre- 
sented in JOI)- Last, but not least, the investigation of the 
exclusive final states at large momenta allows for tests of 
our theoretical understanding of form factors within the 
framework of perturbative QCD. Beyond the intrinsic in- 
terest in this reaction, these studies may provide impor- 
tant clues for the interpretation of exclusive decays of B- 
mesons, a topic of evident importance for the extraction 
of CKM matrix elements. 

Measurements of the cross section for electron-positron 
annihilation were traditionally performed in the scanning 
mode, i.e. by varying the beam energies of the collider. 
The recent advent of and B-meson factories with enor- 
mous luminosities allows us to exploit the radiative re- 
turn to explore the whole energy region from threshold 
up to the nominal energy of the collider. Photon radiation 
from the initial state reduces the cross section by a fac- 
tor 0(a/n). However, this is easily compensated by the 
enormous luminosity of 'factories' and the advantage of 
performing the measurement over a wide range of ener- 
gies in one homogeneous data sample |11| (for an early 
proposal along these lines, see J2). in principle the reac- 
tion e + e~ — > 7 + hadrons receives contributions from both 
initial- and final-state radiation. Only the former is of in- 
terest for the radiative return; the latter has to be elim- 
inated by suitably chosen cuts. The proper analysis thus 
requires the construction of Monte Carlo event generators. 
The event generator EVA was based on a leading order 
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treatment of ISR and FSR, supplemented by an approx- 
imate inclusion of additional collinear radiation based on 
structure functions and included two-pion and four- 
pion final states |13| . Subsequently the event generator 
PHOKHARA was developed; it is based on a complete 
next-to-leading order (NLO) treatment of ISR [T1I15II161 
IT7lli8lll9ll2()| . In its version 2.0 it included ISR at NLO 
and FSR at leading order (LO) for tt + tt~ and fi- 
nal states (and four-pion final states without FSR in the 
formulation described in detail in 

Recent preliminary experimental results indeed demon- 
strate the power of the method and seem to indicate that a 
precision of one per cent or better is within reach [2111221 

a further improvement of our theoretical understanding 
seems to be required. In the present work we therefore 
present a new version (version 3.0) of PHOKHARA, which 
allows for the 'simultaneous' emission of one photon from 
the initial and one photon from the final state, requiring 
only one of them to be hard. This includes in particular the 
radiative return to 7r + 7r _ (7) and thus the measurement of 
the (one-photon) inclusive 7r + 7r~ cross section. The new 
version of PHOKHARA and various tests of its technical 
precision will be presented and several suggestions for a 
model-independent extraction of the form factor will be 
discussed. 

The issue of photon radiation from the final states is 
closely connected to the question of 7r + 7r~ (7) contribu- 
tions to a M (for related discussions, see e.g. [71l3"3"] 1. Soft 
photon emission is clearly described by the point-like pion 
model. Hard photon emission, with E 1 > 0(100 MeV), 
however, might be sensitive to unknown hadronic physics. 
We will therefore study the size of virtual, soft and hard 
corrections separately, and argue that contributions from 
the hard region, above 100 MeV, are small with respect to 
the present experimental and theory-induced uncertainty. 

The outline of the present paper is the following. In 
section 2 we recall the basic aspects of the radiative re- 
turn in leading order, define strategies to separate ISR 
and FSR through cuts, and discuss their markedly differ- 
ent angular distribution and the characteristic feature of 
ISR-FSR interference. In section 3 we present the estimate 
of 7T + 7r~ (7) contributions to a^. Section 4 is devoted to the 
description of the additional two-photon corrections orig- 
inating from simultaneous ISR-FSR emission and their 
implementation in PHOKHARA 3.0. The impact of these 
new terms on various distributions is described in section 
5. A brief summary and our conclusions can be found in 
section 6. 



2 ISR versus FSR and the radiative return at 
leading order 

Let us, in a first step, recall the basic ingredients for a 
proper exploitation of the enormous luminosity at <P- and 
_B-meson factories through the radiative return. In this pa- 
per we will concentrate on the measurement of the n + ir~ 
final state, accompanied by one or two photons. We shall 
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Fig. 1. Leading order contributions to the reaction e + e 
Tr+Tr^ from ISR (a) and FSR (b). 



start with a discussion of the leading process, with only 
one photon, radiated from the electron (positron) or the 
hadronic system: 



e + e 



ir + (p + ) 7T (p ) 7 



(1) 



For the radiative return this corresponds to the Born ap- 
proximation, and higher order corrections will be discussed 
below. The amplitude of interest describes radiation from 
the initial state (ISR) and is depicted in Fig. (per- 
mutations of photon lines will always be omitted). It is 
proportional to the pion form factor, evaluated at Q 2 = 
(p + +p~) 2 . However, radiation from the charged pions 
(FSR) (Fig. [1J)) evidently leads to the same final state 
and must be suppressed and controlled with sufficient pre- 
cision. A variety of methods, which have already been de- 
scribed in |11II16II18| . will be recalled in the following. 

The fully differential cross section describing photon 
emission can be split into three pieces 



da = rfcTTSR + da 



'FSR 



da 



INT 



(2) 



which originate from the squared ISR and FSR amplitudes 
respectively, plus the interference term. They depend on 
two Dalitz variables, which characterize the energies of 
7r + and 7T _ and the photon, and on the three Euler angles 
describing the orientation of the 7r + 7r~7 production plane 
in the centre-of-mass system (cms). The dependence on 
the azimuthal angle around the beam direction is trivial 
in the case of unpolarized beams. 

The interference term, dciNT, is antisymmetric under 
the exchange of 7r + and tt~ or e + and e~. This allows us 
to test the model for FSR and, furthermore, can even lead 
to a measurement of the amplitude Afsr. The usefulness 
of this method has already been emphasized in |11|. 

Let us, for the moment, concentrate on the charge- 
symmetric pieces. The inclusive photon spectrum, sep- 
arated according to ISR and FSR, is shown in Fig. [21 
The discrimination between ISR and FSR is based on 
the markedly different angular distribution of the photon. 
Without any assumption on the nature of FSR, the cross 
section can be written in the form 



Q 



da 



dQ 2 d cos 9 n 



4^ 2 ^ 

6s 4 

(s 2 + Q 4 ) 1_ 

s(s - Q 2 ) 1 - cos 2 1 



Q 



2 - 



2s 



+Ai(s,Q 2 )cos 2 



A 2 {s,Q 2 )sin 2 e i , (3) 



Henryk Czyz et al.: The radiative return at and B-factories: FSR at next-to- leading order 



3 



120 



$ 100 



a so 

60 

t 

^ 40 



fcs 



20 



a 






• * 

• 




1 


1 

.02 GeV 








• 

• 


• 




• 









» 

■ 


• 

• 

■ 

• 




• 






* 

8 


■ 

■ 

■ 


. i Vis 

— ■ . . — *-• 


• 

R+FSR "* " 

***** 
SR only.' 

■■■■■■ 

Ri I SR 

r 1 > 165° 
r =t < 140° 


***** 



s 

t* 

It* ■ 

■■■ 


— -»— ■ 
* . 

■ 

■ 


■ 


■ * 
■. i 

■ 

■-. 

i IS 

& 1 < 15° c 

40° < 
i 



0.2 0.3 0.4 0.5 0.6 0.7 0. 

Q 2 (GeV 2 ) 



8 0.9 



0.25 



0.2 



P4 

^2 0.15 



Pi 

PLh 0.1 



0.05 - 



-0.05 



"i 1 1 r 



"1 1 1 

1.02 GeV 



no cuts 



x 10 



# < 15° or ft, > 165° 
40° < n ± < 140° 
J I I I I I I 



0.5 
<jT^ 0.45 
^ 0.4 
0.35 

a 

0.3 

'b 0.25 

+ 

f 0.2 

j~ 0.15 
b 

•ra 0.1 
0.05 




0.002 



c 






10.5^ 


Z GeV 



































xs 

























m 


ISR+ 


FSR 






> 
















xs 








< X 





















******* 



0.2 0.4 0.6 0.8 
Q 2 (GeV 2 ) 

T 



1.2 



0.0015 - 



-<£ 0.001 - 



Pi 

j£i 0.0005 - 



- -0.0005 - 



0.2 0.3 0.4 



0.5 0.6 0.7 
Q 2 (GeV 2 ) 



0.8 0.9 1 



-0.001 




0.4 0.6 0.1 
Q 2 (GeV 2 ) 



1.2 



Fig. 2. Contributions to the inclusive photon spectrum from ISR compared to ISR plus FSR at y/s = 1.02 GeV (a) (without 
and with cuts) and \fs — 10.52 GeV (c). In (b) and (d) the relative contribution of FSR is plotted (multiplied by a factor of 
10 for y/s = 1.02 GeV, when cuts are applied). 



where Q 2 is the invariant mass of the ir + , ir system, 9 1 
is the photon polar angle and /3(s) = yl — 4m 2 / s. Terms 
~ m 2 1 s are neglected for simplicity, while they are present 
in the program. Fitting the observed photon distribution 
to this form allows an unambiguously separation of ISR 
and FSR. The composition of the angular distribution of 
photons from the two components is displayed in Fig. [3] 
for three photon energies. 

The second option is based on the fact that FSR is 
dominated by photons collinear to tt + or 7r~, ISR by pho- 
tons collinear to the beam direction. This suggests that 
we should consider only events with photons well sepa- 
rated from the charged pions and preferentially close to 
the beam |11II16II18| . The contamination of the events with 
FSR is reduced to less than five per mille (Fig.JSJs). Con- 
sidering just the inclusive photon distribution, without in- 
clusion of any cut would necessarily have led us to wrong 
conclusions. For higher beam energies, the LO FSR is nat- 



urally suppressed and at y 7 ! = 10.52 GeV it is entirely 
negligible (Fig. |2Ji)- 

The interference term is asymmetric under the ex- 
change of 7r + and 7r~ or electron and positron: 

damT(p + ,p~) = -dawT(p~ ,P + )- (4) 

This gives rise to a forward-backward asymmetry of the 
inclusive pion distribution 1 1 1 1 . In kinematical regions where 
FSR is relatively more pronounced, the interference term 
is relatively large (Fig.gJ:, Q 2 between 0.84 and 0.85 GeV 2 ) 
and vice versa (Figs.^Ji and^Js). Note that all these con- 
siderations refer to the low energy analysis of KLOE. At 
high energies, e.g. at B-factories, ISR and FSR are clearly 
separate for the 7r + 7r~7 final state. 

Similar tests were already performed by KLOE |23j. 
where it was demonstrated that inclusive photon asym- 
metries are in excellent agreement with the EVA MC , 
where FSR emission is modelled by point-like pions. Given 
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Fig. 3. Angular distribution of photons with E 1 £ (314, 319) 
MeV (a), E-y G (221, 226) MeV (b) and £ 7 G (93,98) MeV (c), 
separated according to ISR and FSR. The curves represent the 
fits according to Eq. . 



Fig. 4. Contributions of ISR, FSR and INT terms to pion 
polar angle differential cross section for three different ranges 
of Q 2 . 



enough statistics, the charge-asymmetric piece of the cross 
section can even be investigated for every point in the 
Dalitz plane, e.g. as a function of Q 2 and the angle be- 
tween 7r + and 7 separately. Keeping these variables fixed 
and varying the angle between photon and electron, the 
relative amount of the ISR amplitude can be varied dras- 
tically, while the FSR amplitude stays roughly constant. 



Alternatively, we could formulate the angular distri- 
butions in terms of independent helicity amplitudes for 
FSR, which depend on the two Dalitz variables only, and 
deduce information on these amplitudes from their inter- 
ference with the ISR amplitude. In the absence of any 
model for FSR beyond scalar QED (sQED) and in view 
of the fact that sQED will presumably provide a satisfac- 
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tory description at DA$NE energies we shall not dwell 
further on this subject. 

Let us summarise the main points of this 'leading or- 
der' discussion: 

1. The cross sections for photon emission from ISR and 
FSR can be disentangled as a consequence of the marked 
difference in angular distributions between the two pro- 
cesses, Eq. (01 • This observation is completely gen- 
eral and does not rely on any model like sQED for 
FSR. This allows us to measure the cross section for 
er(e + e - — > 7* — > tt + tt~j) directly for fixed s as a func- 
tion of E*y , an important ingredient in the analysis of 
hadronic contributions to a M , as discussed in Section 
Alternatively we can employ suitable angular cuts to 
separate the two components. 

2. Various charge-asymmetric distributions can be used 
for independent tests of the FSR model amplitude, 
with the forward-backward asymmetry as simplest ex- 
ample. Typically the charge asymmetry is large in the 
region where ISR and FSR are comparable in size and 
small when the separation is clean and simple from 
general considerations. A typical case for this second 
possibility is the radiative return at y/s = 10.52 GeV, 
at the i3-meson factories, where the 7r + 7r~7 final state 
is completely dominated by ISR. 



3 Higher order contributions to the 
anomalous magnetic moment of the muon 

In leading order the invariant mass of the 7r + 7r~7 system 
is restricted to the cms energy ^fs. If the relevant angu- 
lar distributions, as predicted by sQED, coincide with the 
experimental analysis, the validity of this model at lower 
energies is highly plausible. Nevertheless, a measurement 
of 7* — > 7r + 7r~7 for variable y/s is desirable for an inde- 
pendent cross check. In SectionFJ]we will demonstrate that 
this is indeed possible with the radiative return, if events 
with two photons in the final state are investigated. This 
is in fact one of the motivations for extending the event 
generator PHOKHARA to events with simultaneous ISR 
and FSR. However, before discussing this issue in detail, 
an investigation of 7* — * 7r + 7r~7, the corresponding vir- 
tual corrections and the relevance of these amplitudes to 
the analysis of the muon anomalous magnetic moment 
is in order. 

Let us concentrate on the one-particle irreducible ha- 
dronic contributions 

2 />00 7 

«r Lo = h / - r w < ( 5 ) 

Stt z J 4m 2 s 

with the familiar kernel K and the R ratio defined through 

i?( S )cx|(0|J M |had,( 7 ))| 2 . (6) 

R(s) can be obtained from the cross section for electron- 
positron annihilation into hadrons, after correcting for ISR 
and for the modifications of the photon propagator, which 



are sometimes (not quite correctly in the time-like region) 
summarized as 'running' of the fine structure constant. 

At the present level of precision of hadronic contribu- 
tions to Qu, roughly half to one per cent, photonic cor- 
rections to final states with hadrons start to become rel- 
evant. Qualitatively the order of magnitude of this effect 
can be estimated either by using the quark model with 
m u k, md ~ m s sa 180 MeV (adopted to describe the 
lowest order contribution [S3]), or with m u ps rrid ~ m s w 
66 MeV (adopted to describe the lowest order contribution 
to a(Mz)) or by using 7r + 7r _ as dominant intermediate 
hadronic state plus photons coupled according to sQED. 
The three estimates (see also [7]) 

<5a M (quark,7,m 9 = 180 MeV) = 1.880 x 10~ 10 , 
<5a M (quark,7,m g = 66 MeV) = 8.577 x 10~ 10 , 

Sa^+ir- ,j) = 4.309 x lO" 10 , (7) 

are comparable in magnitude and begin to be relevant 
at the present level of precision. This order-of-magnitude 
estimate suggests that a more careful analysis is desirable. 

In view of the dominant role of 7r + 7r _ for the evaluation 
of a^ ad,LO we will from now on concentrate on this specific 
final state. In Section 5 we will argue that the measure- 
ment of i?(7r + 7r~(7)) is indeed possible with the technique 
of the radiative return. Before entering this discussion, a 
precise definition of the objects of interest is required. We 
will rely on the smallness of the fine structure constant 
and only consider NLO amplitudes and rates. In leading 
order the contribution to R and a^ ad ' LO is given by the 
square of the pion form factor. In most experiments the 
issue of contributions from events with additional photon 
emission from the hadronic system and the influence of 
these photons on cuts has simply been ignored (for a first 
step towards inclusion and control of these effects from 
the experimental side, see for related theoretical dis- 
cussions, see [7|Eni)- These effects indeed are only relevant 
in next-to-leading order. 

To proceed to 0(a) for point-like particles such as 
muons or electrons, we simply evaluate corrections to the 
vertex function from virtual photon exchange, which must 
be combined with real radiation to arrive at an infrared- 
finite result. For pions we must follow a different strat- 
egy. The evaluation of virtual corrections involves the full 
hadronic dynamics and seems to be difficult, if not im- 
possible, but is also unnecessary. Instead, we separate the 
inclusive rate i?(7r + 7r~(7)) into a part containing the vir- 
tual plus real photonic corrections i? v+R (7r + 71-^(7), E cut ), 
with energies up to a cutoff E cut chosen such that the 
point-like pion model describes real emission, say up to 
50 or even 100 MeV, and a remainder R h (tt + tt^^, E cut ) 
with hard photons only. The separation between the two 
configurations will depend on the cutoff _E cut , which has 
to be chosen small for this dependence to be given by the 
familiar result for point- like pions. 

The strategy for the exclusive measurement now pro- 
ceeds as follows: 

The contribution from i? v+R (7r + 7r~ (7), E cnt ) as a func- 
tion of _E cut can be measured by analysing final states 
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Let us now estimate the contributions from hard pho- 
ton radiation to . In Fig. [S] we display the integrand 
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Fig. 5. Differential contribution to a)i ad ' 7 from 7r + 7r~7 inter- 
mediate states for different cutoff values compared with the 
complete contribution (virtual plus real corrections, labelled 
'inclusive') evaluated in sQED (FSR), as well as with the con- 
tribution from the 7v + tv~ intermediate state. 
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where hard-photon radiation is suppressed by suitable cho- 
sen cuts on energy and collinearity of the pions. The E cut 
dependence of R b -(tt + tt^j, E cut ) can be measured by col- 
lecting 7r + 7r~7 final states (and correcting of course for 
ISR) as discussed in Section [2] Alternatively we may cor- 
rect for (unmeasured) hard photon events by employing a 
model like sQED, which has to be checked experimentally 
— at least for a few selected kinematic configurations. 



-7-o^ ad,7 (J5 c 
ds M 



3tt 2 



-K(s) i? H (s,£ c 



(8) 



for E cut = 10, 100 and 200 MeV as a function of ^/l = 
m(7r + 7r~7) between the threshold and 2 GeV. The result 
is compared with the complete sQED contribution, as de- 
rived from point-like pions, and the lowest order contri- 
bution from 7r + 7r~. The integrated result is displayed in 
Fig.© 

Using sQED as a model, one finds that contributions 
from hard photon radiation, with a cut at 100 MeV, are 
still small with respect to the present uncertainty 5a{j ad,LO = 
7 x 10 -10 [EE]. Let us emphasize that cuts on the pho- 
ton energy around 50 MeV or below might well lead to 
important shifts in a^. Of course we have assumed that 
hard radiation is not grossly underestimated by sQED, an 
assumption to be tested by experiment. As we will demon- 
strate in the next section, such tests are indeed feasible in 
experiments based on the radiative return. 



4 FSR at next-to-leading order and the event 
generator PHOKHARA 

As discussed in Section |21 the radiative return allows us 
to exploit the enormous luminosity of <2>- and B-meson 
factories for the measurement of the hadronic cross section 
over a wide range of energies. On the basis of the leading 
order treatment, we have demonstrated that the analysis 
of final states with one photon allows us to determine the 
pion form factor for arbitrary Q 2 and the cross section for 
7r + 7r~ plus a hard photon from FSR at y/s. 

Intuitively it should also be possible to exploit the ra- 
diative return for the extraction of FSR for arbitrary in- 
variant mass sfs' of the 7r + 7r~7 system through the reac- 
tion 

e + e~ — > 77*(— > 7r + 7r~7) . (9) 

The implementation of this two-step process into PHO- 
KHARA and the question of how to separate the corre- 
sponding amplitude from double photon emission from the 
initial or final state are the main topics of this section. 

Let us first describe the main physics ingredients and 
assumptions of the new version of PHOKHARA (PHO- 
KHARA 3.0). The complete set of diagrams relevant to 
the reaction 

e + e~ — ► 7r + 7r~77 , (10) 

is displayed in Fig.0 We require at least one photon (with 
four momentum k\) to be hard, i.e. with energy E\ > 
E™ ln , where £™ in is typically around 5 to 10 MeV for 
DA$NE energies and significantly larger (several GeV) for 
B-meson factories. We only display typical diagrams and 
omit permutations. However, we explicitly distinguish the 
two cases, where the hard photon is emitted either from 
the electron-positron line or from the hadronic system. 
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(c) (d) 

Fig. 7. Typical amplitudes contributing to the reaction 
e + e~ — > 7r + 7r~77. For two photons emitted either from the 
electron/positron or the hadronic system. Only one represen- 
tative is displayed. 




(c) (f) 

Fig. 8. Typical amplitudes describing virtual corrections to 
the reaction e + e~ — > 7t + 7t _ 7. Permutations are omitted. 

Diagrams relevant to virtual corrections to the reaction 
with one photon, 

e + e~ — > 7T + 7r~7 , (11) 

are shown in Fig. |S1 they contribute through their inter- 
ference with the Born amplitude (Fig.^). At this point no 
assumption about the form of the 7* — > 7r + 7r~7 coupling 
or the structure of the virtual corrections at the hadronic 
vertex is imposed. 

Just as in leading order, charge-symmetric and asym- 
metric terms are present in the differential distribution. 
The latter are particularly sensitive to the FSR ampli- 
tude and may serve for tests of the FSR model or even for 




(F) 



Fig. 9. Charge-even infrared-finite combinations of ampli- 
tudes. The photon labelled T' is hard. Diagrams written in 
reverse order are to be understood as complex conjugate. 



an independent determination of this amplitude for arbi- 
trary s' . In the present upgraded version of the program, 
only charge symmetric terms will be included, which are 
sufficient for the present purpose. 

In Fig. those combinations of (gauge-invariant sets 
of) amplitudes involving real and virtual radiation are dis- 
played, which lead to infrared-finite charge-even combina- 
tions. The first (A) corresponds to the radiative correc- 
tions to the leptonic tensor, and has been derived in 
117) . It is the crucial element of PHOKHARA (versions 
1.0 and 2.0) and can be used for the simulation of arbi- 
trary hadronic final states. The second combination (B) 
corresponds to the two-step process described in Eq. © 
above and is implemented in the present, newest version 
of PHOKHARA (version 3.0). Before discussing the ac- 
tual implementation, let us discuss the remaining terms. 
The next combinations, denoted by (C) and (D), describe 
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radiative corrections to pure FSR, with the leading pro- 
cess depicted in Fig. The hard photon originates from 
the hadronic system, and the radiative corrections (vir- 
tual, soft, and hard radiation) are either attached to the 
hadronic (C) or the leptonic (D) system. As discussed in 
Section [SJ the leading order process can be controlled at 
DA<J>NE through suitable cuts and FSR can be reduced to 
a level of less than 1%. At the energy of 10.52 GeV, i.e. at 
_B-mcson factories, it is entirely irrelevant. Therefore we 
do not include these radiative corrections to an already 
small contribution in the program. 

A similar argument applies to the remaining terms, de- 
noted by (E) and (F). These describe radiative corrections 
(real and virtual photons) to the interference between hard 
ISR and hard FSR. They will be neglected on the basis of 
the same arguments as before. 

In the following, we denote by ISRNLO the contribu- 
tions to the cross section coming from ISR only calculated 
at NLO, i.e. those from Figs.Hk andEK- IFSLO includes 
also the contribution from LO FSR (diagrams in Figs-QJi, 
UJd and|3Y)- Finally, IFSNLO includes on top of that ini- 
tial plus final state emission at NLO (diagrams in Figs.^Ji, 
and |33). Let us emphasize that these are gauge- 
invariant subsets throughout. 

Let us now describe the implementation of the two- 
step process Eq. (J5J, denoted IFSNLO, in detail. Since 
we are only interested in final states with masses of the 
7r+7r~7 system below « 1 GeV, we expect sQED to pro- 
vide a good description of the 7* — > 7r + 7r _ 7 amplitude. 
The virtual corrections are similarly modelled by sQED. 
They are combined with soft radiation, with E2 < -E 2 ut , 
to arrive at an infrared- finite result. The effective form 
factor is, therefore, introduced through the substitution 



2\|2 



\F\Q 2 )\ 



£f ut ) 



\F{Q 2 
\F°(Q 2 )\ 2 



n v+s 



TTTCUt 

• h 2 



)) ,(12) 



with ?? V+S (Q 2 , E 2 ut ) given below. For the generation of 
real radiation the cutoff energy obviously refers to the 
7T + 7r~7 centre-of-mass frame. 

It should be emphasized that only the combination 
\F(Q 2 , iJ™*)! 2 is physically observable and relevant to the 
analysis of and a(mz), of course after adding real ra- 
diation, which can either be measured, or for photon en- 
ergies up to 0(100 MeV), modelled by sQED through the 
reaction 



"(Pi>A e +) + e (>2,A e -) -> 

7r+((7i) + ir~(q 2 ) + 7(fci,Ai) + j(k 2 , A 2 ) 



(13) 



The helicity amplitudes responsible for real radiation 
(Fig. CJd and c) are given by the sum of 12 Feynman am- 



plitudes. Adopting a notation similar to 16 , they read 
M H (\ \ \ \ ) — ( 4 ™) 2 / 

JW IFSNL0l A e+' A e-,^l>^2j — ^ | 

v](pi,X e +) A u I (p 2 ,\ e -) +v} I (pi,X e +) B u n (p2,X e -) 
+ (fci<->fca) , (14) 

where the index H indicates that both photons are hard, 
and 



A = 



(e*(k 1 ,X 1 )-k+ ~2e*(k 1 ,\ 1 )- Pl ) J- 



B 





2ki ■ pi 


J- (2e*(fci, Ai 


)-P2-k+e*{k l ,\ l )-) 




2k\ ■ P2 


(e*(fci,Ai)+*; 


- -2e*(k u X 1 )-p 1 )J+ 




2ki ■ pi 


J+ (2e*(k u Ai 


) ■P2-Ke*{k u \ 1 )+) 


2ki ■ p 2 



(15) 



with 



(16) 



(17) 



Q = Pi + P2 - fci = 9i + 92 + k 2 , s = Q' 
and J M being the current describing the 7r + 7r _ 7 final state 
J» = iF 2 „(s') D" , (18) 

where 



£>" = ( qi + k 2 - q 2 



„ qi ■ e*(fc 2 , A 2 ) 



qi ■ k 2 



^ 92 ■£ (fc2,A 2 ) 

+ (?2 + K2 — q\) : 

92 • h 



2e^(fc 2 ,A 2 ) , (19) 



and k ± 1 are 2x2 matrices, which can be found in |16|. 
We do not include contributions from the radiative decay 
<j) — > 7r + 7r _ 7. Their effect is small, can be well controlled 
and was discussed in 

In the present version of the program, one of the pho- 
tons is assumed to be visible, at least in principle, and 
only the photon emitted from final states is allowed to 
be soft (Fig. |5j3). For photon energies 11;^ < £i 7 < E% 
(or Wyfs < E 2j < El) the amplitude used consists of 6 
diagrams only. E v c is the threshold above which we can 
'observe' the photon — typically 10 MeV for DA$NE 
and 100 MeV for B-factories. If both photons are hard, 
Ei 7 , E 2l > E v c the sum of 12 diagrams is used. 

The corresponding virtual plus soft photon corrections 
can be written as 



r1rr V+S - - r7 V+S (V F CU M Hrr (0 > (VI 



(20) 



where t&jjg^ ^ s ^ ne l ea( hng order e + e — * tt + it 7 cross 
section, with the photon emitted off the initial leptons 
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only, and 



77 V+ V,£: 2 cut ) = -2 



l + ft 2 



2/3 



log(t) + 1 



log(2 W ) 



ft 

l + ft 2 



a 2 (i-ft 2 

log(t)-2 + log' 



2/3 



log(t) log 



/3(l + /3) 



log 
2Li 2 



1 + (1-ft 

^rj l0g 

2/3 



where w = E^/Vs*, and 



2/3 
2Li 2 



1-/3 
2/3 



/? = y/1 - Aml/s' , i 



1-ft 
1 + /3 



(21) 



(22) 



The function ?7 V+S (s', -E 2 ut ) ^ s °f course the familiar 
correction factor derived in [37] for the reaction e + e~ — » 
7r + 7r _ 7 in the framework of sQED (see also |SHj)- For the 
present case, s' corresponds to the squared mass of the 
7r + 7r~7 system, and for virtual and soft photon emission, 
s' ~ Q 2 , the soft photon cutoff is defined in the 7r + 7r~7 
rest frame. 

Correspondingly the cross section for the reaction 



7(fci)7*(-> vr+vr 7(fc 2 )) , 



(23) 



after integration over the angles and energy (from to 
the kinematical limit) of 7(^2) , is given by 



j_H 

acr IFSNLO 



(24) 



with 



r, H (s',E^) 



l + ft 2 



ft 



Li 2 1 



C(2) 



log 2 (t) 



2 

(1 - ft 2 )ft, 



\0g{t m - t) 

2(1 -ft 2 ) 



2/3 3 (l-/%) 



1 - 2 



l + ft 2 



7ft 2 



log( T 



, log(^m) 
+ (1+/? 2 ) 



-2 
log 



1-/3 2 



/? 2 



4(/? 2 -/3 r 2 J 



log(i T 



,(l + /3) 2 (l-/32„)y ' 2 
-21og(l-f m i) , 

where i and ft are denned in Eq. (|2"2"|) . and 



(25) 



/3 m = 4/I- 



Ami 



l-A, 



Qfn 1 "I" A' 

Here is the maximum value of Q 2 



Or. 



2En 



and i? 2 ut is again defined in the 7* rest frame. This for- 
mula will be useful for tests of the Monte Carlo event gen- 
erator discussed below, and £" 2 ut is still arbitrary at this 
point. The same formula can also be used to test FSR at 
LO with s' — > s and £" 2 ut set now in the e + e~ cms frame. 

For E% nt small, w = E% ut /V? < 1, the function rf 1 
reduces to 



V H (s',E^ 
(l + ft 2 ) 



log(2tu) 



1 + /3 2 

/3 



log(t) 



ft 



log(t 2 )log(l-i)+Li 2 (l-i 2 ) 



(l-/3 2 )(3 + /3 2 ) 
4/3 3 



log(t) + 21og 



l-ft 2 
4/3 2 



3 
2^2 



(26) 



Adding virtual, soft (Eq. (J2U) and hard (Eq. ffify) cor- 
rections, the familiar correction factor ET7lI3??f7| 



') = 

- 3 log 
1 



1 + /3 2 

/? 



4Li 5 



log 



I- ft 
l + ft 
l+ft 



2Li 2 



4 V 



41og(/3) 



1 + /3 

/3 2 ) 2 - 2 
3 (1+/3 2 ) 



1-/3 



-21og(/3)log 



1-/3 
1 + /3 
1 



■/3 



log 



1 + /3 
1-/3 



3 log 



1-/3 
1-ft 2 



2/3 2 



(27) 



is recovered. 

Using Eqs. i(2"T)l and (|23|l the implementation of FSR 
in combination with ISR is straightforward. To match 
hard, soft and virtual radiation smoothly, the energy cut- 
off (£ l 2 Ut ) has to be transformed from the rest frame of 
the 7r + 7r~7 (emitted from the final state) system to the 
laboratory frame (e + e~ cms frame) (S™ m ). In fact it is 
necessary to recalculate the soft photon contribution, as 
the cut on Q 2 depends in the latter case on the angle 
between the two emitted photons and now 



r7 V + fa ( S , S ',m 



l + ft 2 



2ft 



log(t) + 1 



log(2w 
l + ft 2 



-7^ ^ 4 
s' — s \ s 



ft 

l + ft 2 
2/3 

1 



log(t) - 2 + log 



1-ft 2 



log 
2Li 2 



- log(t) log 



4 

ft(l + ft) 



2/3 
2ft 



log 



2/3 



l + ft 



2Li 2 - 



1-ft 
2/3 



— 7T 

3 



(28) 



where ft is denned in Eq. (|2~2"j) and w — E™ m /y/s . 

A number of tests were performed to ensure the tech- 
nical precision of the new version of PHOKHARA. The 
square of the matrix element summed over polarizations 
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Table 1. Total cross section (nb) for the process e + e~ — > 
7r + 7r~7(7) for different values of the soft photon cutoff. One of 
the photons is required to have an energy larger than 10 MeV 
(100 MeV) for s/s = 1.02 GeV (10.52 GeV). No further cuts 
applied in columns 2 and 3. In column 4, the pion angles are 
restricted to 40° < 8 W ± < 140° and the photon(s) angles to 
7 < 15° or 9 1 > 165°. 

w Js =1.02 GeV 10.52 GeV 1.02 GeV 



10 
10 
10 



40.992 (5) 
41.013 (6) 
41.018 (7) 



0.1606 (1) 
0.1607 (2) 
0.1607 (2) 



20.988 (1) 
20.993 (1) 
20.995 (2) 



of the final particles and averaged over polarizations of 
the initial particles was calculated with FORM gU| us- 
ing the standard trace method. External gauge invariance 
was checked analytically when using the trace method and 
numerically for the amplitude calculated with the helic- 
ity amplitude method. The two results for the square of 
the matrix element summed over polarizations, were com- 
pared numerically. The code based on the result from the 
trace method was written in quadrupole precision to re- 
duce cancellations. The code based on the result obtained 
with the hclicity amplitude method, uses double-precision 
for real and complex numbers and is now incorporated in 
the code of PHOKHARA 3.0. Agreement of 13 significant 
digits (or better) was found between both codes. The sen- 
sitivity of the integrated cross section — with and without 
angular cuts — to the choice of the cutoff w can be de- 
duced from Table ^ an( l Fig- QDI For simplicity the same 
separation parameter w was chosen for ISR and FSR cor- 
rections. Choosing w = 10 -4 or less, the result becomes 
independent of w. The tests prove that the analytical for- 
mula describing soft photon emission as well as the Monte 
Carlo integration in the soft photon region are well im- 
plemented in the program. Having analytical expressions 
for ry H (s, E™ ln ) (Eq. (|25|l ) we can test also the implemen- 
tation of hard photon emission from the final state. The 
results of the tests are collected in Figs. II II and 1121 where 
FSR at LO obtained from PHOKHARA is compared with 
the analytical result for the corresponding cross section: 



cfsrO) 



mm 

7 



a H/ 

- V (s, 

7T 



77111111 

^1 



) Ms) 

7ra 2 /3 3 (s) 



I F 2l! {s) | 2 , (29) 



for a fixed value of the photon energy cutoff and differ- 
ent values of the e + e~ cms energy (Fig. llljl or fixed value 
of the e + e~ cms energy and various values of the photon 
energy cutoff (Fig. I12[). As one can see the technical pre- 
cision of that part of PHOKHARA is much better than 1 
per mille. The IFSNLO part was also tested against the 
analytical result 



r)(s') = 



d V+S+H 
""IFSNLO 



(s,s')/W 



§ da^(s')/ds' 



(30) 



The results of the tests are shown in Fig.^| The agree- 
ment of r\ calculated by PHOKHARA and the analytical 
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Fig. 10. Comparison of the Q 2 distribution for two values of 
the soft photon cutoff (w = 10" 4 vs. 10~ 5 ) for ^/s = 1.02 GeV 
(a,b) and y^s — 10.52 GeV (c). One of the photons was required 
to have energy > 10 MeV (for ■ s /s = 1.02 GeV) and > 100 MeV 
(for y/s — 10.52 GeV). Angular cuts were applied only on plot 
(b), where pion polar angles are restricted to 40° < 6 7T ± < 140° 
and photon polar angles to (9 7 < 15° and # 7 > 165°. 



result (Eq. J23) is at the level of a few per mille. However, 
since the contribution to the cross section of that part is 
multiplied by a/ir ~ 1/400, the actual technical precision 
of the IFSNLO PHOKHARA cross section is at the level 
of 10~ 5 . 



Henryk Czyz et al.: The radiative return at and B-factories: FSR at next-to- leading order 



11 



0.001 



0.0005 



-0.0005 - 



-0.001 



e^e — y tt^tt 7 

1 1 1 — 1 — 1 1 1 1 J 



a(PHOKHARA) 
a(anal) 



IV 



10 



-4 



FSR at LO 



0.001 



e T e — > 7r 7r 7 



10 



100 



(GeV) 



Fig. 11. Comparison between the LO FSR cross section calcu- 
lated analytically (Eq. (gpj) and calculated by PHOKHARA 
for fixed value of w = E™ in /yfs. 



5 FSR contribution to the hadronic cross 
section and its measurement via the radiative 
return method 

We shall study now the impact of the new corrections 
on various distributions. Before entering the discussion, 
let us recall the meaning of various abbreviations, which 
will be used in the following. ISRNLO corresponds to the 
ISR cross section calculated at NLO without any FSR. 
IFSLO includes in addition FSR at LO. Finally, IFSNLO 
stands for ISR and FSR at NLO as implemented in the 
new version of PHOKHARA (version 3.0). 

Let us start the discussion for a cms energy of 1.02 
GeV, relevant for the KLOE experiment. The IFSNLO 
correction to the cross section from graphs in Fig. |5}3 is 
relatively big at low Q 2 , if no cuts are applied (Fig. 114b'). 
Below the p resonance, they grow from zero at the res- 
onance to 10% of the IFSLO cross section near the pro- 
duction threshold, while they remain small above the p 
resonance. This is due to the fact that the ISR leading 
to the p meson is strongly enhanced. Subsequently the p 
decays into 7r + 7r~7, with a large contribution from the re- 
gion where Q 2 , the invariant mass of the 7r + 7r~, is low. 
Of course this is smeared by the width of the p, but the 
above discussion remains valid and the contribution from 
the newly implemented diagrams, through the reaction 
e + e~ — > 7/)(— > 77r + 7r - ), is sizeable. 

The effect of several 'standard cuts' at y/s = 1.02 GeV 
is shown in Figs. 114b and!14b. The sensitivity of the newly 
implemented contributions to these cuts can be exploited 
to test, control or experimentally eliminate FSR at NLO. 
It is relatively easy to find cuts that lower the correction 
to a level of 2%-3%. In fact all the cuts which were used to 
eliminate FSR at LO are also effective here. The standard 
KLOE cut on the track mass [2H| reduces FSR further, to 
less than 1% for most of the Q 2 range, apart of the high 
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Fig. 12. Comparison between the LO FSR cross section calcu- 
lated analytically (Eq. ^) and calculated by PHOKHARA 
for a fixed value of = 1.02 GeV. 



values of Q 2 , where the corrections remain at the level of 
2%-3% (Fig.Glt, lower curve). 

The behaviour of the additional contribution at 10.52 
GeV is qualitatively similar. Without cuts the relative en- 
hancement amounts to approximately 4% (Fig. 115b ) , and 
this remains unchanged after cuts on pion and photon 
angles, which correspond to the detector acceptance, are 
introduced (Fig. H5b). As stated before, a large fraction 
of the excess is due to the feed-down from the p-meson in 
the two step process e + e _ — > "/p(-^ TT + n~"f). Assuming, 
that hard photons, say above 100 MeV, can be identi- 
fied, we try to identify these events by assigning photons 
with small angles relative to n + or ir~ (say 15°) to the 
hadronic system when M p — r p < -J (p + + p_ + p 7 ) 2 < 
M p + r p , thus moving these events from Q 2 — (p + +p-) 2 
to Q 2 = (p + + p- +Pj) 2 . For a sample generated with 
PHOKHARA 3.0 without FSR at NLO the distribution 
remains practically the same (Fig. 115b ). while including 
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FSR at NLO the excess in the low Q 2 region is signifi- 
cantly reduced (Fig. I15t0. In fact, part of the remaining 
enhancement (about 1 to 2% for Q 2 between 0.2 and 0.1 
GeV 2 ) just corresponds to the effect of multiplication 
with the factor —rj(s'). 

The FSR corrections originating from the 'two-step' 
process and their effect on cuts necessarily introduce some 
model dependence into the extraction of the pion form 
factor. The detection and analysis of the complete final 
state — charged pions as well as photons — is of course 
optimal for the study of this phenomena. However, the 
present KLOE analysis is based on the measurement of 
the 7r + and n~ momenta only. Even this partial informa- 
tion allows to construct distributions, which are sensitive 
to final states with tt + n~ and two photons. As an alter- 
native to the straightforward study of the Q 2 distribution 
and its Monte Carlo simulation it is thus even possible to 
test the model governing the production of 7r + 7r~ in con- 
junction with two photons and constrain or even measure 
the process 7* — > 7r + 7r~(7), an important ingredient for 
the precise prediction of a^ ad,L °. 

The response of dcriFSNLo/^Q 2 to different cuts on 
M 2 , the invariant mass of the recoiling photonic system 
(Fig. I16|) will be an important observable. Contributions 
with double emission from the initial state and from mixed 
ISR-FSR will be affected by this cut. However, the cut de- 
pendence is markedly different for duiFSNLO /dQ 2 (which 
includes both contributions) and dcrisRNLo/^Q 2 (which 
includes configurations from ISR only) as shown in Fig. 1171 
The confirmation of the M 2 dependence as predicted by 
PHOKHARA, as well as the verification of charge asym- 
metric distributions (see Section 2) would provide addi- 
tional support for the ansatz of emission from point-like 
pions or allow to test alternative models. 



6 Conclusions 

Measurements of the pion form factor through the radia- 
tive return offer the unique possibility for improved pre- 
dictions for the muon magnetic moment and the electro- 
magnetic coupling a(Mz)- Contributions to the dispersion 
integral from intermediate states with ir + ir~ and a pho- 
ton start to become relevant at the present level of pre- 
cision and can be measured in reactions leading to 7r + 7r~ 
in conjunction with two photons. An upgraded version of 
PHOKHARA (version 3.0), which includes simultaneous 
emission of photons from the initial and the final state, 
has been presented. The two step process e + e~ — > 7p(— > 
7r + 7r _ 7) leads to a notable enhancement of events with 
low mass of the 7r + 7r - system. Various cuts are described, 
which allow to control this effect, to identify these events, 
correct the distribution and test the model for final state 
radiation. 
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